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The contribution of self-induced charge currents of metamaterial media to photon effective rest 
mass is discussed in detail in the present paper. We concern ourselves with two kinds of photon ef- 
fective rest mass, i.e., the frequency- dependent and frequency-independent effective rest mass. Based 
on these two definitions, we calculate the photon effective rest mass in the left-handed medium and 
the 2TDLM media, the latter of which is described by the so-called two time derivative Lorentz 
material (2TDLM) model. Additionally, we concentrate primarily on the torque, which is caused 
by the interaction between self-induced charge currents in dilute plasma {e.g., the secondary cosmic 
rays) and interstellar magnetic fields (ambient cosmic magnetic vector potentials), acting on the 
torsion balance of the rotating torsion balance experiment. 



I. PHOTON EFFECTIVE REST MASS DUE TO SELF-INDUCED CHARGE CURRENT 

In is well known that in some simple electromagnetic media such as electron plasma, superconducting media and 
Lorentz dispersive materials (inside which the time-harmonic electromagnetic wave is propagating), the self- induced 
charge current density J [1] is proportional to the magnetic vector potential A. Thus the interaction term /xqJ • A in 

2 2 

the Lagrangian density is therefore transformed into an effective mass term — ^ "^""^ of electromagnetic fields [2], 
which is analogous to the rest mass term in the Londons' electromagnetics for superconductivity, Ginzberg-Landau 
superconductivity theory and Higgs mechanism. Note that in these media, the magnetic permeability^ n = 1, so it 
is easy for us to obtain J (and hence the effective rest mass of photon) by solving the equation of motion of charged 
particles acted upon by electromagnetic waves. However, for some artificial composite electromagnetic metamaterials, 
which are described by, e.g., the two time derivative Lorentz material (2TDLM) model (in which both the electric 
permittivity and the magnetic permeability are of frequency dependence and therefore of complicated form) and some 
uniaxially (or biaxially) anisotropic media with permittivity and permeability being tensors, the above formulation 
may be not helpful in obtaining the effective rest mass of electromagnetic fields. In these cases, however, we have 
shown that the following formula [2] 

^ = (l-n>^ (1.1) 

which arises from the Einstein-de Broglie relation, can be applicable to the above problem. In the expression (1.1), 
n, h and c denote the optical refractive index, Planck's constant and speed of light in a vacuum, respectively. 

The two time derivative Lorentz material (2TDLM) model first suggested by Ziolkowski and Auzanneau [3,4] is 
a generalization of the standard Lorentz material model. Ziolkowski has shown that this type of 2TDLM medium 
can be designed so that it allows communication signals to propagate in the medium at speeds exceeding the speed 
of light in a vacuum without violating causality [5]. For these artificial metamaterials, the model encompasses the 
permittivity and permeability material responses experimentally obtained [6]. The x-directed polarization and y- 
directed magnetization fields in such materials would have the following forms [5] 



*This paper is concerned with the self-induced charge currents and some related physically interesting topics such as photon 
effective rest mass due to medium dispersion and its effects in torsion balance experiments. It will be submitted nowhere else 
for the publication, just uploaded at the e-print archives. 

^Electronic address: jqshen@coer.zju.edu.cn 

^Note that here for the magnetic properties in superconductors, we adopt the current viewpoint rather than the magnetic- 
charge viewpoint (where the permeability can be viewed as jj. — 0). 
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My + T^My + UjlMy = Lolx^ny + LO^X'^^^Hy + Xj^Hy, (1.2) 



where Xq™; X^'™ ^ind xl^'™ represent, respectively, the couphng of the electric (magnetic) field and its first and second 
time derivatives to the local electric (magnetic) dipole moments. cUp, r°, and ojq can be viewed as the plasma 
frequency, damping frequency and resonance frequency of the electric (magnetic) dipole oscillators, respectively. Thus 
the frequency-domain electric and magnetic susceptibilities are given 

^ {^-)\l+ia:u;;x}-u;'x; ^ (O^xg' + i^c.-xg'-a^^xy 3 

So, the refractive index squared in the 2TDLM model reads 

(u;) = 1 + x'^ {u^) + X" (c^) + X° (c^) x"^ (^) • (1-4) 
Thus it follows from (1.1) that the photon effective rest mass squared in the 2TDLM model is of the form 

^eff = ^ [X^ (^) + X"" {^) + X" ('^) X'^ (u^)] ■ (1-5) 

As an illustrative example, we apply Eq.(1.5) to the case of left-handed media, which is a kind of artificial compos- 
ite metamaterials with negative refractive index in the microwave frcqiunic-y region (10^ Hz). More recently, this 
medium captured considerable attention in various fields such as condensed matter physics, materials science, applied 
electromagnetics and optics [6-11]. The most characteristic features of left-handed media are: (i) both the electric 
permittivity e and the magnetic permeability /i arc negative; (ii) the Poynting vector and wave vector of electromag- 
netic wave propagating inside it would be antiparallel, i.e., the wave vector k, the electric field E and the magnetic 
field H form a left-handed system; (iii) a number of peculiar electromagnetic and optical properties such as the rever- 
sal of both the Doppler shift and the Cherenkov radiation, anomalous refraction, amplification of evanescent waves 
[12], unusual photon tunneling [13], modified spontaneous emission rates and even reversals of radiation pressure to 
radiation tension [7,8] arise. All these dramatically different propagation characteristics stem from the sign change of 
the group velocity. In what follows we calculate the photon effective rest mass squared by making use of Eq.(1.5) and 
the expressions 

2 2 

e{u,) = l- M(a>) = l- , ■ r (1-6) 

for the effective dielectric parameter and magnetic permeability, where the plasma frequency ojp and the magnetic 
resonance frequency ujq are in the GHz range {e.g., LOp = 10.0 GHz, = 4.0 GHz [14]); 7 and F represent the 
damping parameters. The parameter F may often be chosen < F < 1, for example, F = 0.56 [14]. By using 
{lo) = e (lo) ji (w), we obtain [15] 



"<"' = '- J(i+trj + \^--^+-\ ^ (") 

It follows that the mass squared m^fj is obtained 
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+ (1-^ + -) +Fu;^(l-^ + - 



(l.E 



Note that for the electron-plasma medium, where the electromagnetic parameters 7 = F = 0, F = 0, one can arrive 
at the following familiar formula 

mis = (1.9) 

with the plasma frequency squared uj^ = (here e, rUe and N respectively denote the electron charge, mass and 
electron number density in this electron plasma). 
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It should be emphasized that there is a difference in the definition of photon effective rest m.ass between two types 
of experimental schemes of testing photon mass^. In the discussion of wave dispersion and potential variations of the 
speed of light with frequency, which are based on the Maxwellian equations {rather than the Hamiltonian or Lagrangian 
density of interacting electrorn,agnetic system), the effective rest mass squared of photons is expressed by Eq.(1.9). 
Historically, this problem has been discussed by many investigators, e.g., Feinberg who considered the possibility of 
the variation of the speed of light with frequency in the sharply defined optical and radio pulses from pulsars and took 
into account in more detail the observed variation of arrival time with frequency for the radio waves attributed to 
the interaction with interstellar electrons [17]. However, for those experimental schemes of testing photon rest mass, 
which are based on the Hamiltonian density {rather than the Maxwellian equations), the effective rest mass of photons 
due to media dispersion (i.e., the interaction of wave with charged particles in environmental dilute plasma such as 
the secondary cosmic rays) is two times that expressed in Eq.(1.9). The torsion balance experiment is just this type 
of schemes, which will be considered in Sec.V. 

However, it should be noted that in Lorentz dispersive materials, m^^ in (1.8) depends strongly on uj, and sometimes 
possesses an imaginary part (due to the damping parameter and resonance frequency) and the physical meanings of 
miff is therefore not very apparently seen. In what follows we will define a physically meaningful effective rest mass 
of photons, which is independent of the frequency cj, and will explain a method by which the frequency-independent 
effective rest mass of photons in the 2TDLM media is extracted from n'^ {oj). Note that the above non-physical 
disadvantage of frequency-dependent mass in (1.8) will be avoided in the frequency-independent effective rest mass. 

II. EXTRACTING FREQUENCY-INDEPENDENT EFFECTIVE MASS FROM REFRACTIVE INDEX 

SQUARED 

In this framework of definition of frequency-independent effective rest mass, the photon should be regarded as acted 
upon by a hypothetical force field, namely, it behaves like a massive de Broglie particle moving in a force field. It 
is verified in the following that, for the case of de Broglie particle in a force field, one can also consider its "optical 
refractive index" n. According to the Einstein- de Broglie relation, the dispersion relation of the de Broglie particle 
with rest mass mo in the potential field V (x) agrees with 



{oj - cj>f 



2„4 



(2.1) 



with 



It follows that 



CO 



1 - 



-2^ + 



(2.2) 



which yields 
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Compared with the dispersion relation ^ = , one can arrive at 



(2.3) 



1 



2„,2 



(2.4) 



which is the square of "optical refractive index" of the de Broglie particle in the presence of a potential field </>. 

If a photon is permitted to possess a frequency-independent effective rest mass mefi, then the expression (2.4) for 
the case of photons in media is rewritten 



^These two types of experimental schemes are as follows: one is based on the Hamiltonian density of electrodynamics and the 
other based on Maxwellian equations. 
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n^iu,) = l-^-2^+(^) . (2.5) 

In what follows, an approach to frequency- independent effective mass extracted from the optical refractive index 
squared r? (w) is presented. It is assumed that r? (w) can be rewritten as the following series expansions 



-'H=E^. (2-6) 
fe 

then one can arrive at r? (w) = ^0^^ namely, 

(w) = ... + + a-2 + ct-io; + oqw^ + a+iw^ + ... . (2.7) 
w 

Compared (2.7) with (2.5), it is apparent that the frequency-independent effective rest mass squared of the 
photon can be read off from v? (w), namely, the constant term on the right-handed side on (2.7) is related only to 
m^ff, i.e., 

a-2 - — ^2~- y^-^> 

Generally speaking, the optical refractive index (uj) is often of complicated form, particularly for the metama- 
terials. So, we cannot extract frequency-independent m^g- immediately from (uj). In the following, this problem is 
resolved by calculating the limit value of (w) w^, where u! tends to oo or zero. We first consider the large frequency 
behavior of {u))w^. It follows from (2.7) that {w)!^"^ at high frequencies {i.e., u — > oo) is of the form 

(a;)a;^— >a_2 + fl-iw + ao'^^ + tt+iw^ -|- a+2W* -h ... . (2.9) 

This, therefore, means that (w) at high frequencies {u — > oo) can be split into the following two terms 

(w) (constant term) -|- (divergent term). (2.10) 

So, we can obtain rn^g from the following formula 

2 4 

lim Tn^ (w) - (divergent term)] = - . (2.11) 

Note that the constant term in (w) can also be extracted by considering the low-frequency (zero-frequency) 
behavior of (w). It is certain that in this case the obtained constant term is the same as (2.11). Apparently 
{u})ui^ at low-frequencies {i.e., w — > 0) behaves like 

2/ \ 2 , C'-5 , a_4 a_3 n^ 

n {ivjuj ^...-\ -\ -\ \-a-2, (2.12) 

where the sum term (... -|- -|- + ^^) is divergent if w approaches zero. Thus, subtraction of the divergent term 
from {uj) lv^ yields 

2 4 

lirn^ [n^ {(jj) - (divergent term)] = - "^^^^ . (2.13) 

In the next section, we will calculate the frequency-independent effective rest mass of photons in the 2TDLM model 
by making use of the method presented above. 



III. FREQUENCY-INDEPENDENT EFFECTIVE REST MASS OF PHOTONS IN THE 2TDLM MODEL 

According to the definition of the 2TDLM model, the frequency-domain electric and magnetic susceptibilities are 
expressed by Eq.(1.3) and the refractive index squared is given in Eq.(1.4) 

In order to obtain the frequency-independent m^fj from r? {iS), according to the formulation in the previous section, 
the high-frequency {lo oo) behavior of x'' (w) o;^ should first be taken into consideration and the result is 
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(3.1) 



Note that —iuLUpX^ + A7 is the divergent term of x** (t^) at large frequencies. Subtracting the divergent term 
from x** (t^) t^^, we consider again the high-frequency behavior {u 00) and obtain another divergent term ioj^VXj: 



Thus, after subtracting all the divergent terms from x** {oj)u)'^, we obtain 



;rx^-(a;0)'x; + (rTx; 



(3.2) 



(3.3) 



jirn^ [x' (t^) - (-iww^x^ + w^x;) - iw^r^x;] = - (t^;j xs 

In the same fashion, we obtain 

Jim [x'" (a;) - {-iuJu^^X^ + u;\^) - iu^^T^x"^] = - X^ - ^F^x^ - W)' X^T + (r'")' X^] (3-4) 

for the magnetic susceptibility x™ {'^)- In what follows we continue to consider the term x'' (t^) X™ (i^) 011 the right- 
handed side of (1.4). In the similar manner, the infinite-frequency limit of x** {<^) X™ (^) is given as follows 



hm x^(w)x"^(cc)a>' 



-X, 



[iolfxl - u^lT^xl {<f X"^ + in' X", 



K)'x; 



^ Xfi - K^o ) Xj +{^ ) X-i 



(3.5) 



Hence, insertion of (3.3)-(3.5) into (2.11) yields 



(1 + x:^) H)'xi - '^^r^x^ - Hf x; + in' x; 



+ (1 + x;) «)' x^ - <r-x^ - x-^ + {rn' x". 



(3.6) 



Thus the w-independent effective rest mass squared m^g is extracted from {ui). Apparently, m^jj in Eq. (3.6) 

2 4 

does not depend on oj. The remainder on the right-handed side of Eq. (2.5) is (w) — 1 -|- ™|^^2 ; which is equal to 
it) ~ ' wiisre (j) is defined to be ^. For simplicity, we set 
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A(a;)=n^(a;)-1 + ^, 



(3.7) 



and then from the equation 



UJ 



(3.8) 



one can readily obtain the expression for the hypothetical potential field 4> (w) is (j) (w) =ijO 1 — -\/l -|- A (w) 



, or 



y (w) = /iw 1 - Vi + A (w) 



(3.9) 



Many metamatcrials which have complicated expressions for electric and magnetic susceptibilities such as (1.3) are 
generally rarely seen in nature and often arises from the artificial manufactures. More recently, one of the artificial 
composite metamatcrials, the left-handed medium which has a frequency band (GHz) where the electric permittivity 
(e) and the magnetic permeability [^) are simultaneously negative, has focused attention of many authors both 
experimentally and theoretically [6,9-11,15,19]. In the left-handed medium, most phenomena as the Doppler effect. 
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Vavilov-Chcrcnkov radiation and even Snell's law are inverted. In 1964^, Veselago first considered many peculiar 
optical and electromagnetic properties, phenomena and effects in this medium and referred to such materials as 
left-handed media [7], since in this case the propagation vector k, electric field E and magnetic field H of light 
wave propagating inside it form a left-handed system. It follows from MaxwcUian curl equations that such media 
having negative simultaneously negative e and n exhibit a negative index of refraction, i.e., n = —y/eji. Since 
negative refractive index occurs only rarely, this medium attracts attention of many physicists in various fields. In 
experiments, the negative e and /i can be respectively realized by using a network (array) of thin (long) metal wires 
[19] and a periodic arrangement of split ring resonators [9]. A combination of the two structures yields a left-handed 
medium. Compared the permittivity and permeability (1.6) of left-handed medium with the electric and magnetic 
susceptibilities (1.3), one can arrive at 

xl = 1, = ^'p: x)3 = X7 = 0, Lo^o = 0, r'' = -7, 

X^ = -F, xa=xf = 0, <=a;o, r- = -r. (3.10) 

So, it is readily verified with the help of (3.6) that the w-independent effective rest mass of photons inside left-handed 
media is written as follows 

^ = (l-F)a;^ + F(a;o^-r^). (3.11) 

It follows that Eq.(3.11) is a restriction on the electromagnetic parameters F,cOp,uio,J^ and F in the electric per- 
mittivity and magnetic permeability (1.6), since m^jj > 0. Insertion of the experimentally chosen values of the 
electromagnetic parameters in the literature into Eq. (3.11) shows that this restriction condition is satisfied. For 
instance, in Ruppin's work [14], F, Wp, ljq, 7, F and F were chosen F = 0.56, Up = lO.OGHz, uq = 4.0GHz, 7 = 0.03a;p, 
r = 0.03a;o5 which indicates that the inequality m^g > is satisfied. 



IV. SOME APPLICATIONS OF FREQUENCY-INDEPENDENT EFFECTIVE REST MASS OF 

PHOTONS 

In this section we consider the applications of the concept of frequency-independent effective rest mass of photons 
to some optical and electromagnetic phenomena and effects. 

(i) As was stated above, since the frequency-independent effective rest mass of photons in electromagnetic media 
is related close to the coupling parameters of the electric (magnetic) fields to the electric (magnetic) dipole systems, 
and the damping and resonance frequency of the electric (magnetic) dipole oscillators, it contains the information 
on the interaction between light and materials. It is shown that in certain frequency ranges this effective rest mass 
governs the propagation behavior of electromagnetic wave in media, and by using this concept one can therefore 
consider the wave propagation somewhat conveniently. For example, in the left-handed media whose permittivity and 

permeability are expressed by (1.6), the potential field V{(j) approaches —Ti^F{ujQ — Wp), which is constant, as the 

frequency tends to zero. This, therefore, means that at low frequencies the light propagation is governed mainly by 
the frequency-independent effective rest mass of photons. 

Likewise, in the case of the potential field V{ijj) being vanishing where the frequency squared of electromagnetic 
wave is 

02 = L , (4.1) 

the photon can be regarded as a "free" massive particle with the rest mass being expressed in (3.11). In the frequency 
region near Q., sine the photon is a "quasi-free" particle, the light propagation behavior and properties is rather simple 



^Note that, in the literature, some authors [20,21] mentioned the wrong year when Veselago suggested the left-handed media. 
They claimed that Veselago proposed or introduced the concept of left-handed media in 1968 or 1964. On the contrary, the 
true history is as follows: Veselago's excellent paper was first published in Russian in July, 1967 [Usp. Fiz. Nauk 92, 517-526 
(1967)]. This original paper was translated into English by W.H. Furry and published again in 1968 in the journal of Sov. 
Phys. Usp. [7]. Unfortunately, Furry stated erroneously in his English translation that the original version of Veselago' work 
was first published in 1964. 
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compared with other frequency ranges in which V{uj) is nonvanishing or large. For the left-handed media, the relevant 
frequency parameters wq and iVp of which is about GHZ, the frequency f2 in Eq.(4.1) that leads to V{Q) = is just 
in the frequency band (GHZ) where the permittivity and permeability are simultaneously negative. It follows that 
when the effects and phenomena associated with negative index of refraction occur in left-handed media, the wave 
propagation {e.g., scattering properties) and dispersion properties of the light are not as complicated as those in other 
frequency ranges, since in the negative- index frequency band, the optical refractive index squared simulates nearly 
the plasma behavior, particularly in the region near Q. 

(ii) Historically, possibility of a light pulse with speed greater than that (c) in a vacuum has been extensively 
investigated by many authors [5,22-27]. Ziolkowski has studied the superluminal pulse propagation and consequent 
superluminal information exchange in the 2TDLM media, and demonstrated that they do not violate the principle of 
causality [5]. He showed that in the 2TDLM model, both the phase and group velocities, i.e., 

lim v(u>) ~ — r (4-2) 

of light at high frequencies will exceed the speed (c) of light in a vacuum, so long as — 2 < < 0, where x-y denotes 
the coupling coefficients of the second time derivatives of electric (magnetic) fields to the local electric (magnetic) 
dipole motions [5]. 

Although in the 2TDLM model the medium will exhibit a possibility of superluminal speeds of wave propagation, 

the photon velocity (rather than group velocity of light) is by no means larger than c. This may be explained as 
follows: according to Eq.(3.9), the kinetic energy of a photon is £ = hui^l + A(w). If the light frequency cu tends 
to infinity, then it follows from Eq.(3.7) that £ — > nUw. Since the photon momentum p = the velocity v of the 
photon with a frequency-independent rest mass approaches c, i.e., the particle (photon) velocity v = ^ ^ c. Thus 
we show that in the artificial composite metamaterials such as 2TDLM media the photon velocity does not exceed 
the speed of light in a vacuum, no matter whether the phase (and group) velocity of light inside it is a superluminal 
speed or not. 

(iii) The frequency-independent rest mass of photons is useful in discussing the plus and minus signs of phase shifts 
in light scattering inside media. The phase shifts and its signs contain the information on the scattering between 
wave and potentials. According to Eq.(3.9), if A{uj) is positive (negative), then the interaction energy V{uj) between 
light and media is negative (positive). For this reason, it is possible for us to determine the signs (plus or minus) of 
phase shifts in light scattering process by calculating A(a;) and the frequency-independent rest mass meg. Now with 
the development in design and fabrication of left-handed materials, the zero-index (e ^ 0,/i ^ 0) materials receive 
attention in materials science and applied electromagnetism [28]. For this kind of media with the optical refractive 
index n ~ 0, from Eq.(3.7) it follows that 

^('^) - - 1' (4-3) 

2 4 

i.e., A(a;) is approximately equal to ™|"2 — 1 and we can determine the sign by comparing the light frequency with 

2 

the photon effective rest mass moff. This, therefore, means that if the frequency u) is larger than HhM^^ then A < 
and consequently the phase shift of light wave due to scattering acquires a minus sign; and if the light frequency oj is 

2 

less than HMfc_ ^ \\^q^ A > and consequently the phase shift acquires a plus sign. 



V. EFFECTS OF SELF-INDUCED CHARGE CURRENTS IN THE TORSION BALANCE EXPERIMENT 

In this section, we consider the potential effects of self-induced charge currents arising in the torsion balance 
experiments [29,30]. The Lagrangian density of electrodynamics reads 

= - ^^7^,*^/^ + 1^0 Jn\, (5.1) 

where — (^^)^ with m^, h and c being the photon rest mass, Planck's constant and speed of light in a free 

vacuum, respectively, and jiQ denotes the magnetic permeability in a vacuum and the sum over the repeated indices 
is implied. It follows from (5.1) that the canonical momentum density is written 

dC 

Tr{x) = --E{x). (5.2) 
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In the electron plasma, due to the conservation law of the canonical momentum density, i.e., ^(rricV 



have V — — - 

me 

density J = A^ev 



C with C being a constant velocity. By using the formula Ti 
- —^^A + NeC, one can arrive at (in SI) 



E A 



eA) = 0, we 
£ and electric current 



n 



where Jq 



1 

p and Aq 



^+B^ + M^A^ + 2^A^. 



1 



me 



V- E- 



Jo 



+ JqAo - NeC ■ A, 



(5.3) 



are respectively the electric charge density and electric scalar potential, and eo represents 
the electric permittivity in a vacuum and c = J-^ . Note that here ( V • E — 4^ ) results from both the mass 



2 a I Jo 



term -^Aq and the Gauss's law V • E = -/x^Aq + ^. 

It follows from (5.3) that the total effective rest mass squared of electromagnetic fields in electron plasma is given 
as follows 



2 2 , cl^oNe^ 
/"tot = < + 2 



(5.4) 



It is worthwhile to point out that according to the Ampere-Maxwell-Proca equation V x B = /ioeo^E + hqJ — n^A 
and the consequent VxB = ^qCo^E- [/x^ + '^"J^f ] A + /xq iVeC , the /x^t should be |J,^^^^ = M7 + '^"J^f rather than that 
in (5.4). This minor difference between these two /x^ot results from the derivative procedure applied to the Lagrangian 
density by using the Euler-Lagrange equation. Although it seems from the Ampcre-Maxwell-Proca equation that 
Mtot = + ^°m^ ' factor-2 in (5.4) cannot be ignored in calculating the torque acting on the torsion balance 
due to the photon effective rest mass, which will be confirmed in what follows. In the classical electromagnetics, the 
following familiar formulae are given 

V X B = /xqJ, rud = nr'^I, / = J • S, r = ma x B 



and 



VxA = B, ad = TTr^ — , $ = B-S 
Mo 



(5.5) 



(5.6) 



where V x B = /xqJ holds when the electric field strength E, electric current density J are vanishing in the Ampere- 
Maxwell equation. In Eq.(5.5), = nr'^I means that a current loop of radius r carrying current / gives rise to a 
magnetic dipole moment m^; S in / = J • S denotes the area vector through which the electric current density vector 
J penetrates; r = nid x B means that an electric-current loop with a magnetic dipole moment mrf, immersed in a 
magnetic field B, experiences a torque r = x B. Eqs.(5.6) means in the torsion balance experiment, a toroid coil 
contains a loop of magnetic flux $ which acts as a dipole source of magnetic vector potential via V x A = B with 
the magnetic field B within the toroid as the source term. 

It follows from (5.3) that in the Hamiltonian density Ti. of electromagnetic system, there are |B^ and 



y : ~ „j . A^ (i.e., i (/Xtot-A-) ), which provide US with some useful insights into calculating the torque produced 
by the interaction of magnetic dipole vector potential moment with the ambient cosmic magnetic vector potential 
acting upon the toroid in the torsion balance experiment. For convenience, the Equations (5.6) may be rewritten as 
follows 



V X (/xtotA) = /xo 



^ MtotB 



Mo 



Mtotctd = wr 



f /xtot£\ Mtot'I' _ / MtotB 
\ Mo / ' Mo \ Mo 



•S. 



(5.7) 



In order to calculate the torque acting on the torsion balance, we compare the equations and expressions in (5.7) 
with those in (5.5) as follows: 





the familiar case 


the case due to photon mass 


inH 


IB^ 


i (^totA)"^ 


Equations 


V X B = /xqJ 


Vx(MtotA)=Aio(^) 


dipole moment 


rrid = nr^I 


Mtota. = 7rr-2(^) 


flux 


/ = J • S 


A»tot* _ / MtotB \ g 

Ho \ 1^0 J 


torque 


T = nid X B 


f =? 



*Table 1 
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Thus it follows from the analogies in the above table that, ? = (/itotad) x (/itot-A.), i.e., the torque f = ^i^^^S-d x A. 
Here the torque r arising from the interaction between the dipole field of magnetic potentials (resulting from the 
electric current / in the toroid windings) with the ambient vector potential A acts on the torsion balance. It should 
be noted again that here /i^^t = + 2^^^^ rather than fj,^^^ = fj3, + ^'^^f-- 

Note that there exists a term /xoiVeC in the Ampere-Maxwell-Proca equation V x B = ^qco^E — [^^ + ]A + 
HqNcC Docs this term has influences on the torsion balance? In what follows wc will discuss this problem. 

In the similar fashion, we consider the familiar Hamiltonian density — M • B, which describes the interaction of 
the magnetic moments with the magnetic field B. According to the Ampere-Maxwell-Proca equation with /xoeo^E 

and A being ignored, we have V x B = ncNeC, which can be rewritten as V x B = (/Uo/xtot) ("^fv)' 

meanwhile, H = —NeC ■ A can also be rewritten asH = — ("^f^^ ■ (MtotA). 

It is known that in the sufficiently large space filled with homogeneous magnetic moments with volume density 
(magnetization) being M, the interior magnetic field strength B produced by homogeneously distributed magnetic 
moments is just /XqM^, i.e., V x A = hqM. or V x (/XtotA) = (nol^tot) M. 

Thus the comparisons between "i/ie familiar case" and the case due to the constant C are illustrated in the following 
tabk\ 





the familiar case 


the case due to constant C 


in n 


IB^ 


i (^totA)"^ 


interaction Hamiltonian density 


n = -MB 


^ = -(tr)-(/^-A) 


Equations 


V X (MtotA) = (AtoA*tot)M 


VxB = (MoMtot)(ff^) 


torque 


f = M x B 


f=? 



*Table 2 



It follows that ? = y^jj^j x (MtotA), i.e., f = NeC x A, which dose not act on the toroid in the torsion balance 
experiment. 

Since we have discussed the related preliminary preparation for the effects of self-induced charge currents, we think 
the order of magnitude of its effects in experiments {e.g., the extra torque due to the photon effective mass acting on 
the toroid in the torsion balance experiment) deserves cvahiation. It is shown by our evaluation that the muons and 
alpha-particles in secondary cosmic rays will contribute an effective rest mass about 10~^^ Kg to the photon, which 
is compared to the newly obtained upper limit on photon rest mass in Luo's rotating torsion balance experiment. 



VI. DISCUSSION AND REMARKS: SELF-INDUCED CHARGE CURRENTS IN DILUTE PLASMA 

Here we discuss the effective rest mass of photon resulting from the self-induced charge currents in the torsion 
balance experiments [29,30]. The self- induced charge current arises mainly from the two sources around us: (i) 
muon (/i) component and alpha-particles in cosmic rays; (ii) decay daughter (alpha-particles) of radioactive radon gas 
(II^Rn, ll'^Rn) in room environment. 

It is known that at the sea level, the current density of muon component in secondary cosmic rays is about 
1 X 10~^cm~^ • s~^ [33]. Assuming that the muon velocity approaches speed of light, the volume density of muon can 



''This may be considered in two ways: (i) if the direction of M is assumed to be parallel to the third component of Cartesian 
coordinate system, then the potentially non-vanishing components of B and H are also the third ones. It follows from the 
Ampere's law V x H = Jfree with Jfree = that H3 is constant, and H3 can be viewed as zero. So, according to the expression 
B = /L(o(H -h M), we obtain B = /itoM; (ii) In accordance with Ampere's circuital law, the magnetic induction B inside a long 
solenoid with N circular loops carrying a current / in the region remote from the ends is axial, uniform and equal to /io times 
the number of Ampere-turns per meter nl, i.e., B — ^onl, where n — with I being the solenoid length. Let 5* denote the 
cross-section area of this long solenoid. The volume density of magnetic moments IS resulting from the carried current / is 
M = ii^lni with the solenoid volume V = SI. This, therefore, means that the volume density of magnetic moments just equals 
the number of Ampere-turns per meter nl, i.e., M = nl. So, we have B = fj,oM. 
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be derived and the result is A^^ = 0.3 x 10 ^ m ^. So, according to Eq.(5.4) {i.e., nies = f y ^^^)) electromagnetic 

wave with wavelength^ A ^ 100 m at the sea level acquires an effective rest mass about 0.3 x 10"^"^ Kg. 

It is believed that the air molecules in the vacuum chamber of torsion balance experiment cannot be ionized by the 
alpha-particles of cosmic rays, because the mean free path of alpha-particle moving at about lO^'m • s~^ in the dilute 
air with the pressure being only 10~^ Pa [29] is too large (more than 10^ m)^. So, the air medium has nothing to 
contribute to the effective rest mass of photons. 

Radon (geRn) possesses the a-radioactivity. Generally speaking, the Rn density in room environment is 10~^m~^ or 
so, which may be the same order of magnitude as particle density in secondary cosmic rays^, and the decay daughter 
(alpha-particles) will therefore give rise to a photon effective mass of about 10"^** Kg. 

Thus the effective rest mass indiiced by muons and alpha-particles in secondary cosmic rays and radon gas can be 
compared to the new upper limit (1.2 x lO"^** Kg) obtained by Luo et al. in their recent rotating torsion balance 
experiment [29] . So, as far as the contribution of ions of cosmic rays to the photon rest mass is concerned, this upper 
limit (1.2 X 10"^'* Kg) becomes just a critical value. We think Luo's experimentally obtained upper limit is of much 
interest and therefore deserves further experimental investigation to get a more precise upper limit on photon mass. 



VII. CONCLUSION 

In this paper, 

i) it is verified that the expression '"'"'^ = (1 — n^)a;^ can be applied to calculating the effective rest mass of 
photons in the artificial composite electromagnetic metamaterials with complicated permittivity and permeability. 

ii) we suggest an approach to obtaining the frequency-independent effective rest mass of photons in certain artificial 
composite metamaterials such as 2TDLM media which possess the complicated electric permittivity and magnetic 
permeability. 

iii) it is shown that the extra torque 2 ^°^^ x A by the self- induced charge currents will also acts upon the 
toroid in the torsion balance experiment whose apparatus is immersed in the dilute plasma such as the cosmic rays 
and radon gas. Here the photon effective rest mass squared is nig = 2 ^°^^ rather than /x^g = . 

The following difference between these two fi^g should be pointed out: the theoretical mechanism of torsion balance 
experiment lies in the Lagrangian or Hamiltonian (5.3) of interacting electromagnetic system, so, the effective rest 
mass squared of electromagnetic fields is fi^g = "2 ^°^^ . However, for those experimental schemes [17,31,32] of photon 
rest mass stemming from the Ampcrc-Maxwell-Proca equation V x B = /ioeo^E— [/i^ + ^"^'^ ]A + ^qNcC, the effec- 
tive rest mass squared of electromagnetic fields is /x^g = rather than /x^g = 2 ^°^'^ . The so-called experimental 
realizations stemming from the Ampere-Maxwell-Proca equation are as follows: pulsar test of a variation of the speed 
of light with frequency [17]; geomagnetic limit on photon mass based on the analysis of satellite measurements of the 
Earth's field [31]; experimental test of Ampere's law at low temperature [32], etc.. In all these experimental schemes, 
the effective mass squared is /x^^ — ^"^^ . 



®Since the mean distance of two muons in secondary cosmic rays at sea level is about 100 m, Lorentz's mean-field formulation 
(Lorentz's dispersion theory, Lorentz's electron theory, 1895) is not applicable to the electromagnetic wave with wavelength 
A <C 100 m. So, the above effective mass formula is no longer valid for this case {i.e., X <^ 100 m). In fact, the electromagnetic 
wave witii wavelength A <C 100 m dose not acquire this effective rest mass. But the ambient cosmic magnetic vector potentials 
(interstellar magnetic fields) with low (or zero) frequencies will truly acquire this effective rest mass. 

^According to some handbooks, the free paths of alpha-particles in air (with pressure 10^ Pa) are 2.5, 3.5, 4.6, 5.9, 7.4, 8.9, 
10.6 cm, corresponding respectively to the energy 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 Mev. So, in the vacuum chamber with pressure 
10~^ Pa in Luo's experiment, the free path of alpha-particles (with energy much more than 10 Mev, even perhaps more than 
100 Mev) in secondary cosmic rays is more than 10*^ times several cm, which means that the air molecules in the low-pressure 
vacuum chamber cannot be easily ionized by the alpharparticles of cosmic rays. 

'^The effective doses of cosmic rays and Rn humans suffer are respectively 0.4 mSv (typical range: 0.3-1.0 mSv) and 1.2 mSv 
(typical range: 0.2-10 mSv). 
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iv) the self-induced charge currents in dilute plasma presented above (i.e., muons and alpha-particles in secondary 
cosmic rays and radon gas) will give arise to an effective rest mass of about lO"^** Kg, which can be compared to the 
new upper limit on photon rest mass obtained in Luo's rotating torsion balance experiment. 
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